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We expect a detectable correlation between two seemingly unrelated quantities: the four point 
function of the cosmic microwave background (CMB) and the amplitude of flux decrements in quasar 
(QSO) spectra. The amplitude of CMB convergence in a given direction measures the projected 
surface density of matter. Measurements of QSO flux decrements trace the small-scale distribution of 
gas along a given line-of-sight. While the cross-correlation between these two measurements is small 
for a single line-of-sight, upcoming large surveys should enable its detection. This paper presents 
analytical estimates for the signal to noise (S/N) for measurements of the cross-correlation between 
the flux decrement and the convergence, (STk), and for measurements of the cross-correlation 
between the variance in flux decrement and the convergence, {(5Tfn). For the ongoing BOSS 
(SDSS III) and Planck surveys, we estimate an S/N of 30 and 9.6 for these two correlations. For the 
proposed BigBOSS and ACTPOL surveys, we estimate an S/N of 130 and 50 respectively. Since 
((8T) 2 k) oc erf, the amplitude of these cross-correlations can potentially be used to measure the 
amplitude of as at z ~ 2 to 2.5% with BOSS and Planck and even better with future data sets. 
These measurements have the potential to test alternative theories for dark energy and to constrain 
the mass of the neutrino. The large potential signal estimated in our analytical calculations motivate 
tests with non-linear hydrodynamical simulations and analyses of upcoming data sets. 
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I. INTRODUCTION 

The confluence of high resolution Cosmic Microwave 
Background (CMB) experiments and large-scale spectro- 
scopic surveys in the near future is expected to sharpen 
our view of the Universe. Arcminute scale CMB experi- 
ments such as Planck [l|, the Atacama Cosmology Tele- 
scope @, Q, the South Pole Telescope 0, i], QUIET 
@ and PolarBeaR Q, will chart out the small scale 
anisotropies in the CMB. This will shed new light on the 
primordial physics of inflation, as well as the astrophysics 
of the low redshift Universe through the signatures of the 
interactions of the CMB photons with large scale struc- 
ture. Spectroscopic surveys like BOSS @, Q and Big- 
BOSS [T(| will trace the large scale structure of neutral 
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gas, probing the distribution and dynamics of matter in 
the Universe. While these two datasets will be rich on 
their own, they will also complement and constrain each 
other. An interesting avenue for using the two datasets 
would be to utilize the fact that the arcminute-scale sec- 
ondary anisotropics in the CMB are signatures of the 
same large scale structure that is traced by the spectro- 
scopic surveys, and study them in cross-correlation with 
each other. In this paper, we present the analytic es- 
timates for one such cross correlation candidate - that 
between the gravitational lensing of the CMB and the 
flux fluctuations in the Lyman-a forest. 

The gravitational lensing of the CMB, or CMB lensing 
in short, is caused by the deflection of the CMBpho- 
tons by the large scale structure potentials [see hll for 
a review]. On large scales, WMAP measurements im- 
ply that the primoridal CMB is well described as an 
isotropic Gaussian random field [Hj]. On small scales, 
lensing breaks this isotropy and introduces a specific form 
of non-Gaussianity. These properties of the lensed CMB 
sky can be used to construct estimators of the deflection 
field that lensed the CMB. Therefore, CMB lensing pro- 
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vides us with a way of reconstructing a line-of-sight (los) 
projected density field from zero redshift to the last scat- 
tering surface, with a broad geometrical weighting kernel 
that gets most of its contribution from the z = 1 — 4 
range [HI, [HI, HH • While CMB lensing is mainly sensitive 
to the geometry and large scale projected density fluc- 
tuations, the Lyman-a forest, the absorption in quasar 
(QSO) spectra caused by intervening neutral hydrogen 
in the intergalactic medium, primarily traces the small- 
scale distribution of gas (and hence, also matter) along 
the line of sight. 

A cross-correlation between these two effects gives us 
a unique way to study how small scale fluctuations in 
the density field evolve on top of large scale over and 
under-densities, and how gas traces the underlying dark 
matter. This signal is therefore a useful tool to test to 
what extent the fluctuations in the Lyman-a flux relate 
to the underlying dark matter. Once that relationship is 
understood, it can also become a powerful probe of the 
growth of structure on a wide range of scales. Since both 
massive neutrinos and dark energy alter the growth rate 
of structure at z ~ 2, these measurements can probe their 
effects. This new cross-correlation signal, should also be 
compared with other existing cross-correlations between 
CMB and LSS that have already been observed and that 
are sensitive to different redshift regimes [H, [l?], EH, EH, 

In this work, we build an analytic framework based on 
simplifying assumptions to estimate the cross-correlation 
of the first two moments of the Lyman-a flux fluctuation 
with the weak lensing convergence k, obtained from CMB 
lensing reconstruction, measured along the same line of 
sight. The finite resolution of the spectrogram limits the 
range of parallel fc-modes probed by the absorption spec- 
tra and the finite resolution of the CMB experiments lim- 
its the range of perpendicular fc-modes probed by the 
convergence measurements. These two effect break the 
spherical symmetry of the fc-space integration. However, 
we show that by resorting to a power series expansion 
it is still possible to obtain computationally efficient ex- 
pressions for the evaluation of the signal. 

We then investigate the detectability of the signal in 
upcoming CMB and LSS surveys, and the extent to which 
such a signal can be used as a probe of neutrino masses 
and early dark energy scenarios. A highlight of our re- 
sults is that the estimated cross-correlation signal seems 
to have significant sensitivity to the normalization of the 
matter power spectrum eg. Consistency with CMB mea- 
surements - linking power spectrum normalization and 
the sum of the neutrino masses - allows to use this cross- 
correlation to put additional constrain on the latter. 

The structure of the paper is as follows. In Sec- 
tion [n] we introduce the two physical observables, the 
Lyman-a flux and the CMB convergence pi A[) , the cross- 
correlation estimators pi B[) and their variances pi C[) . 
Our main result is presented in section III Dl where the 
signal-to-noise ratios are computed. Sect ion Hi El contains 
a spectral analysis of the observables that aims at finding 



the Lyman-a wavenumbers that contribute most to such 
a signal. We focus on two cosmologically relevant appli- 
cations in sections MI Al and IIIIB1 for massive neutrinos 
and early dark energy models, respectively. We conclude 
with a discussion in section HVl 



II. ANALYTICAL RESULTS 
A. Physical Observables 

Fluctuations in the Lyman-a flux 

Using the fluctuating Gunn-Peterson approximation 
[2l| , the transmitted flux T along a los n is related to the 
density fluctuations of the intergalactic medium (IGM) 
<5igm by 



J-(n, z) — exp —A(l 
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where A and f3 are two functions relating the flux fluctu- 
ation to the dark matter overdensities. These two func- 
tions depend on the redshift considered: A is of order 
unity and is related to the mean flux level, baryon frac- 
tion, IGM temperature, cosmological parameters and the 
photoionization rate of hydrogen. A good approximation 
for its redshift dependence is A(z) ~ 0.0023 (I + z) 3 65 
(see [IH). f3 on the other hand depends on to the so- 
called IGM temperature-density relation and in particu- 
lar on the power-law index of this relation (e.g. [23l. |24T|) 
and should be less dependent on redshift (unless temper- 
ature fluctuations due for example to reionization play a 
role, see [HI). For the calculation of signal/noise in the 
paper, we neglect the evolution of A and (3 with redshift. 
While the value of the correlators considered will depend 
on A and /3, their signal-to-noise (S/N) ratio will not. 

On scales larger than about I /i _1 Mpc (comoving), 
which is about the Jeans length at z = 3, the relative 
fluctuations in the Lyman-a flux 5T = (J- — T)j T are 
proportional to the fluctuations in the IGM density field 
IH, H3, HI, EH & We assume that the IGM traces the 
dark matter on large scales, 



-Af35 lGM {n, x) « -A(3S(h, X ). (2) 



The (variance of the) flux fluctuation in the redshift range 
covered by the Lyman-a spectrum is then proportional 
to (the variance of) the fluctuations in dark matter 
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d X (-AP) r S r (n, X ), 
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where the range of comoving distances probed by the 
Lyman-a spectrum extends from Xi to XQ- The r = 1 
case corresponds to the fluctuations in the flux and the 
r = 2 case corresponds to their variance. We stress that 
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the above approximation is valid in linear theory neglect- 
ing not only the non-linearities produced by gravitational 
collapse but also those introduced by the definition of 
the flux and those produced by the thermal broaden- 
ing and peculiar velocities. Note that while the assump- 
tion of "tracing" between gas and dark matter distri- 
bution above the Jeans length is expected in the stan- 
dard linear perturbation theory [3~fl ] , the one between the 
flux and the matter has been verified a-postcriori using 
semi-analytical methods ([HI, and numerical simu- 
lations ([I?], [H, Hil) that successfully reproduce most 
of the observed Lyman-a properties. Furthcmore, non- 
gravitational processes such as temperature and/or ultra- 
violet fluctuations in the IGM should alter the Lyman-a 
forest flux power and correlations in a distinct way as 
compared to the gravitational instability process and to 
linear evolution (e.g. [35l.[36l.[37t). 

Cosmic Microwave Background convergence field 

The effective weak lensing convergence re(n) measured 
along a los in the direction n is proportional to the dark 
matter overdensity 8 through 

K{n,XF) = 2 / dxW L (x,XF) , 1 , 4 

where the integral along the los extends up to a comov- 
ing distance Xf and where W L (x, Xf) = x{Xf~x)/Xf is 
the lensing window function. In what follows we consider 
the cross-correlation of Lyman-a spectra with the conver- 
gence field measured from the CMB, as in Vallinotto et 
al. [38|, in which case xf is the comoving distance to the 
last scattering surface. Note however that it is straight- 
forward to extend the present treatment to consider the 
cross-correlation of the Lyman-a flux fluctuations with 
convergence maps constructed from other data sets, like 
optical galaxy surveys. 

It is necessary to stress here that Eq. (fl) above depends 
on the density fluctuations in the IGM, which in principle 
are distinct from the ones in the dark matter, whereas 
k depends on the dark matter overdensities S. If the 
IGM and dark matter overdensity fields were completely 
independent, the cross-correlation between them would 
inevitably yield zero. If however the fluctuations in the 
IGM and in the dark matter are related to one another, 
then cross-correlating re and <5 J- will yield a non-zero re- 
sult. The measurement of these cross-correlations tests 
whether the IGM is tracing the underlying dark matter 
field and quantifies the bias between flux and matter. 

B. The Correlators 

Physical Interpretation 

The two correlators (STn) and (8!F 2 k) have substan- 
tially different physical meaning: re is proportional to the 



over(under)density integrated along the los and is domi- 
nated by long wavelength modes with k ~ I0~ 2 h Mpc -1 . 
Intuitively re therefore measures whether a specific los is 
probing an overall over(under)dense region. If the IGM 
traces the dark matter field, then by Eq. ([3]) 8 J 7 is ex- 
pected to measure the dark matter overdensity along the 
same los extending over the redshift range Az spanned 
by the QSO spectrum. This implies that 

• (SJ-k) quantifies whether and how much the over- 
densities traced by the Lyman-a flux contribute to 
the overall overdensity measured all the way to the 
last scattering surface. Because both re and 8J- are 
proportional to 8, it is reasonable to expect that 
this correlator will be dominated by modes with 
wavelengths of the order of hundreds of comoving 
Mpc. As such, this correlator may be difficult to 
measure as it may be more sensitive to the calibra- 
tion of the Lyman-a forest continuum. 

• (SJ- 2 k) measures the relationship between long 
wavelength modes in the density and the ampli- 
tude of the variance of the flux. The variance on 
small scales and the amplitude of fluctuations on 
large-scales are not coupled in linear theory. How- 
ever, in non-linear gravitational theory regions of 
higher mean density have higher matter fluctua- 
tions. These lead to higher amplitude fluctuations 
in flux [39|]. Since (SJ- 2 n) is sensitive to this in- 
terplay between long and short wavelength modes, 
this correlator is much more sensitive than (STn) 
to the structure growth rate. Furthermore, because 
ST 2 is sensitive to short wavelengths, this signal is 
dominated by modes with shorter wavelength than 
the ones dominating (SJ-k). As such, this signal 
should be less sensitive to the fitting of the contin- 
uum of the Lyman-a forest. 

Tree level approximation 

In what follows we focus on obtaining analytic ex- 
pressions for the correlations between the (variance of 
the) flux fluctuations in the Lyman-a spectrum and the 
CMB convergence re measured along the same los. From 
Eqs. (J3l [4j) above it is straightforward to obtain the gen- 
eral expression for the signal 

<«-(»)«(»)> = -^y o d Xc a{xc) 

rXQ 

x / d Xq {-A(3) r (8 r (n, Xg )8(n,Xc)). 

(5) 

Since the QSOs used to measure the Lyman-a forest lie 
at z > 2, it is reasonable to expect that non-linearities 
induced by gravitational collapse will not have a large 
impact on the final results. In the following we therefore 
calculate the r = I and r = 2 correlators at tree-level 
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in cosmological perturbation theory. While beyond the 
scope of the current calculation, we could include the 
effects of non-linearities induced by gravitational collapse 
by a pply ing the Hyperextended Perturbation Theory of 
Ref. |5| to the terms in Eq. (©. 

At tree level in perturbation theory the redshift de- 
pendence of the matter power spectrum factorizes into 
p (k,Xc,Xq) = p L(k)D(xc)D(x q ), where P L (k) denotes 
the zero-redshift linear power spectrum and D(x) the 
growth factor at comoving distance x- Furthermore, the 
correlator appearing in the integrand of Eq. ([5]) depends 
on the separation A% = Xq ~ Xc between the two points 
running on the los and in general it will be significantly 
non-zero only when |A%| < A^o ~ 150 /i _1 Mpc. Also, at 
tree level in perturbation theory these correlators carry 
2r factors of D. 1 Using the approximation 



D( X c) 
W l (XcXf) 
a{Xc) 



D( Xg - A X ) 
W L (Xq ~ Ax,xf) 
a(Xg - Ax) 



D(x g ), 

W L (Xq,XF) 

a(Xq), 



(6) 
(7) 
(8) 



we can then write (5 r (h, x q ) S(n, Xc)) ~ (r(^x) D (Xq) 
and trade the double integration (over Xc and Xq) f° r 
the product of two single integrations over A% and Xq- 
Equation ([5]) factorizes into 



-A0) 

-Axo 



3PlQfl r 

2c 2 

dAxU&x) 



Xi a \Xq) 



(9) 



This is the expression used to evaluate the signal. The 
determination of an expression for £ r and of an efficient 
way for evaluating it is the focus of the rest of the section. 



Window Functions 

The experiments that measure the convergence and the 
flux fluctuations have finite resolutions. We approximate 
the effective window functions of these experiments by 
analytically tractable Gaussian function. 

These two window functions act differently: the finite 
resolution of the CMB convergence measurements limits 
the accessible range of modes perpendicular to the los, 
k±, and the finite resolution of the Lyman-a spectrum 
limits the range of accessible modes fc|| parallel to the 
los. This separation of the modes into the ones parallel 
and perpendicular to the los is intrinsically dictated by 
the nature of the obscrvables and it cannot be avoided 
once the finite resolution of the various observational 



1 Notice in fact that even though in the r = 2 case it would be 
reasonable to expect three factors of D, the first non-zero con- 
tribution to the three-point function carries four factors of D 
because the gaussian term vanishes exactly. 



campaigns is taken into account. Because of this sym- 
metry, the calculation is most transparent in cylindrical 
coordinates: k = k»n + fcj_. 

The high-fc (short wavelength) cutoff scales for the 
CMB and Lyman-a modes arc denoted by kc and fc^ 
respectively. Furthermore, we also add a low-fc (long 
wavelength) cutoff for the Lyman-a forest, to take into 
account the fact that wavelengths longer than the spec- 
trum will appear in the spectrum itself as a background. 
We denote this low-fc cutoff by fe;. After defining the 
auxiliary quantities 



P 

A: 2 



h 2 h 2 



k\ h 



2kf 



(10) 
(11) 



the window functions acting on the Lyman-a and on the 
CMB modes, denoted respectively by W a and W K , are 
defined through 



W a (k\\,k L ,ki) 



W K (k ± ,k c ) 



1 



-(fcii/fei) 2 



e -(fcn/^) _ e -(V fc ) 



(12) 
(13) 



where the direction dependence of the two window func- 
tions has been made explicit. 

We determine the values of the cutoff scales as follows. 
For the Lyman-a forest, we consider the limitations im- 
posed by the spectrograph, adopting the two cutoff scales 
fcj, and ki according to the observational specifications. 
For the reconstruction of the CMB convergence map we 
compute the minimum variance lensing reconstruction 
noise following Hu and Okamoto [4l[. We then identify 
the multipole l c , where the signal power spectrum equals 
the noise power spectrum for the reconstructed deflection 
field (for / > l c the noise is higher than the signal). Fi- 
nally, we translate the angular cutoff l c into a 3-D Fourier 
mode kc at the relevant redshift so to keep only modes 
with k < kc in the calculation. Note that if we had used 
the shape of the noise curve instead of this Gaussian 
cutoff, we would have effectively retained more Fourier 
modes, thereby increasing the signal. However, to keep 
the calculations simple and conservative we use the above 
Gaussian window. In what follows, we will present results 
for convergence map reconstructions from the datasets of 
two CMB experiments: Planck and an hypothetical CMB 
polarization experiment based on a proposed new camera 
for the Atacama Cosmology Telescope (ACTPOL). For 
the former, we adopt the sensitivity values of the 9 fre- 
quency channels from the Blue Book [42| ■ For the latter 
we assume a hypothetical polarization based CMB ex- 
periment with a 3 arcmin beam and 800 detectors, each 
having a noise-equivalent-temperature (NET) of 300 fiK- 
yfs over 8000 sq. deg., with an integration time of 3 x 10 7 
seconds. We further assume that both experiments will 
completely cover the 8000 sq. deg footprint of BOSS. 
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FIG. 1: Absolute value of the correlator (SJ-k) along a single line-of-sight as a function of the source redshift z and of the 
length of the measured spectrum Az, for convergence maps recontructed from Planck (left panel, kc = 0.021 /iMpc -1 ) and 
ACTPOL (right panel, kc = 0.064 /tMpc -1 ). The value of the resolution of the QSO spectrum is the one predicted for SDSS- 
III, kh = 4.8/iMpc -1 . To make the physics of structure formation apparent, we turn off the IGM physics by setting A = j3 = 1 
(it is straightforward to rescale the values of the correlator to reflect different values of A and 



Auxiliary Functions 

Because the calculation has cylindrical rather than 
spherical symmetry, the evaluation of the correlators of 
Eq. is more complicated, particularly for r > 1. As 
shown in the appendix, it is possible to step around this 
complication and to obtain results that are computation- 



ally efficient with the adoption of a few auxiliary func- 
tions that allow the integrations in /c-space to be carried 
out in two steps, first integrating on the modes perpen- 
dicular to the los, and subsequently on the ones parallel 
to the los. The perturbative results for the correlators 
are expressed as combinations of the following auxiliary 
functions: 



H m (k\\;kc) = 
L m (k\\;kc) = 



kdk P L {k) 
2tt ml 



k 2 - k 2 
k 2 c 



dk P L {k) ( k2 - k l 
2irk ml 



dk 



f^>(Ax;k c ,k L )= I J 



exp 



exp 



k c 
k 2 -kf 
k c 



k 2 

4 + *h a x 



dk (k 
2tt 1 s 



-k 2 jk 2 



f m (k\\;kc) 
fo(k;oo) 



(14) 



(15) 



with f = {L,H}, (16) 



with /= {L,H}. (17) 



Equations (fl"4|) and (fTS"]) above represent an intermediate 
step, where the integration on the modes perpendicular 
to the los is carried out. Equations (fin)) and (fT?]) are 
then used to carry out the remaining integration over 
the modes that are parallel to the los. 

The symmetry properties of the auxiliary functions are 
as follows. The functions f m are real and even in kti 
regardless of the actual value of m. This in turn implies 
that fm** are real and even (imaginary and odd) in A% 



when n is even (odd). Furthermore, the coefficients /q 
are real and non-zero only if n is even, thus ensuring that 
t; r (Ax) is always real- valued. 

The (8J-k) correlator 

In the r = 1 case it is straightforward to identify 
£i(Ax) with a two point correlation function measured 
along the los. However, the intrinsic geometry of the 
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FIG. 2: Absolute value of the correlator {5T 2 k) along a single line-of-sight as a function of the source redshift z and of the 
length of the measured spectrum Az, for convergence maps recontructed from Planck (left panel, kc = 0.021 /iMpc -1 ) and 
ACTPOL (right panel, kc = 0.064 feMpc^ 1 ). The value of the resolution of the QSO spectrum is the one predicted for 
SDSS-III, kh = 4.8 ftMpc -1 . As before, we set A = /3 = 1 to make the physics of structure formation apparent. 



problem and the inclusion of the window functions leads 
to evaluate this correlation function in a way that is dif- 
ferent from the usual case, where the spherical symmetry 
in fc-space can be exploited. In the present case we have 

6 (A X ) = (A X ; ko,k L )- (A X ; k c , k). (18) 

It is then straightforward to plug Eq. (jT5J) into Eq. © 
to obtain (8T{n)K(h)) ? In Fig. Q]we show the absolute 
value of the cross-correlation of the convergence re of the 
CMB with the Lyman-a flux fluctuations 8T observed for 
a quasar located at redshift z and whose spectrum spans 
a range of redshift Az. The cosmological model used (and 
assumed throughout this work) is a flat ACDM universe 
with f2 m = 0.25, h = 0.72 and as = 0.84 consistent 
with the WMAP-5 cosmology [13] • The left and right 
panel show the results for the resolution of Planck and 
of the proposed ACTPOL experiment. We artificially set 
A = (3 — 1, effectively "turning off" the physics of IGM: 
this choice is not dictated by any physical argument but 
from the fact that it makes apparent the dynamics of 
structure formation. 

The behavior of (5 J- re) shown in Fig. [T] makes phys- 
ical sense. Recall that this correlator is sensitive to 



2 We checked that in the limit where k^ — > oo, kc — > co and fc; — > 
the usual two point correlation function is recovered. Whereas 
one would naively expect that letting k^ = kc and fc ; = would 
lead to recover the usual two point function calculated exploiting 
spherical symmetry in fc-space with a cutoff scale equal to the 
common k^, this is actually not the case. The reason for this is 
that the volume of k- space over which the integration is carried 
out is different for the two choices of coordinate systems. In 
particular, the spherical case always includes fewer modes than 
the cylindrical one. The two results therefore coincide only in 
the fc^ — > oo limit. 



the overdensity integrated along the redshift interval Az 
(spanned by the QSO spectrum) that contributes to the 
CMB convergence. It then increases almost linearly with 
the length of the QSO spectrum Az. It also increases 
if the resolution of the CMB experiment kc is increased. 
An increased value of Az corresponds to a longer Lyman- 
a spectrum, carrying a larger amount of information and 
thus leading to a larger correlation. Similarly, an in- 
creased value of kc corresponds to a higher resolution of 
the reconstructed convergence map and therefore more 
modes - and information - being included in the corre- 
lation. Deepening the source's redshift (while keeping A 
and (3 fixed) on the other hand results in a decrease in 
{8!F re). This fact is related to the growth of structure: 
the spectrum of a higher redshift QSO is probing regions 
where structure is less clumpy and therefore the abso- 
lute value of the correlation is smaller. Finally, once the 
redshift dependence of A is turned on (/3 is only mildly 
redshift dependent) the above result change, leading to 
a final signal that is increasing with redshift. 

We stress here that values of the correlators will be dif- 
ferent when A and (3 are different from unity. Ultimately 
these values should be recovered from a full non-linear 
study based on large scale-high resolution hydrodynami- 
cal simulations. However, numerical studies based on hy- 
drodynamical simulations have shown convincingly that 
for both the flux power spectrum (2-pt function) and flux 
bispectrum (3-pt function) the shape is very similar to 
the matter power and bispectrum, while the amplitude 
is usually matched for values of A and (3 that are dif- 
ferent from linear predictions (see discussion in (43j). In 
this framework, non-linear hydrodynamical simulations 
should at the end provide the "effective" values for A 
and (3 that will match the observed correlators and our 
results can be recasted in terms of these new parameters 
in a straightforward way. 
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The (8J- 2 k) correlator 

The r = 2 case, where the variance of the flux fluc- 
tuation S T 2 integrated along the los is cross-correlated 
with k, is more involved. Looking back at Eqs. (0 
[9]) it is possible to realize that the cumulant correlator 
(S 2 (n, Xq) 8(n, Xc)) = £2 (Ax) corresponds to a collapsed 
three-point correlation function, as two of the 5's refer 
to the same physical point. The evaluation of £2 is com- 
plicated by the introduction of the window functions W a 



and W K . For sake of clarity, we report here only the final 
results at tree level in cosmological perturbation theory, 
relegating the lengthy derivation to the appendix. Let- 
ting 

6(A X )H^c)l.2+ 2(^ 2 ,3, (19) 

and using the auxiliary functions defined in Eqs. (| 14111 7p 
above, it is possible to obtain the following series solution 



(<5 2 <5}i, 2 



+ 



2 E 7 [H^(^x,x q ;kc,k L )-H^(A x , Xq ;k c rk) 

\k L H^\A x , Xq ;k c ,k L ) -kH^(A X ,Xg;kc,k)] [k L L^(A X ,Xg;kc,k L ) ~ kL^(A X ,x q ;kc,k) 
h£\A x , Xq ;k c ,k L )~ h£\a x , Xq ,k c , k)] \l$(A x , Xq ;k c ,k L )- L$(A X , Xq ,k c , k) 



k 2 L L&>{A X , Xq ;k c ,k L )-k 2 L$ (A X , Xq ; k c , k) 



7 

4m 2 



m(2m — 1) 4 
■ kn 



kh L$(Ax,x«;fcc,&L) ~ kL^{Ax,Xq)kc,k) 



7 



L$(A X , Xq ;k Cl k L )~ L<® (A X , x q ,k c , k) 



(20) 



(S 2 q S c ) 2 , 3 = 2 £ 



(-l) m 2 r 



Q -H^\k L )H^\A X ;k c ,k L ) + \k 2 L L^ +1 \k L ) H^ +1 \A x -,k c ,k L ) 



\k 2 L H^ +1 \k L ) 4™ +1) (Ax; k c , k L ) + hi Lt +2 \k L ) 4 m+2) (Ax; k c ,k L ) 

^H^\k L )L^ +2 \A x -k c ,k L ) - ^Lt +2 \k L )Ht\A X ;kc,k L ) + {k L - k) 



(21) 



r 



In Fig. [2] we show the result obtained using the tree 
level expression for (5J- 2 k), Eqs. (|19j|21j) . As before, we 
focus on the physics of structure formation and we turn 
off the IGM physics by setting A = (3 = 1. First, it 
is necessary to keep in mind that (5T 2 k) is sensitive to 
the interplay of long and short wavelength modes and 
it probes the enhanced growth of short wavelength over- 
densities that lie in an environment characterized by long 
wavelength overdensitics. The behavior of (SJ- 2 k) with 
respect to z and Az is similar to that of (5Tk): it in- 
creases if Az is increased or if the QSO redshift is de- 
creased. However, the effect of the growth of structure is 
in this case stronger than in the previous case. This does 
not come as a surprise, as the growth of structure acts co- 
herently in two ways on (5J- 2 n) . Since in a ACDM model 
all modes grow at the same rate, a lower redshift for the 
source QSO implies larger overdensities on large scales 
which in turn enhance even further the growth of over- 
densities on small scales. Thus by lowering the source's 
redshift two factor play together to enhance the signal: 
first the fact that long and short wavelength modes have 



both grown independently, and second the fact that being 
coupled larger long-wavclcngth modes boost the growth 
of short wavelength modes by a larger amount. This de- 
pendence is also made explicit in Eq. © , where we note 
that {8T 2 n) depends on four powers of the growth fac- 
tor. Finally, as before, the higher the resolution of the 
CMB experiment the larger is d(8T 2 k) / dAz. This too 
makes physical sense, as a larger resolution leads to more 
modes contributing to the signal and therefore to a larger 
cross-correlation. 



C. Variance of correlators 

To assess whether the correlations between fluctuations 
in the flux and convergence are detectable we need to es- 
timate the signal-to-noise ratio, which in turn requires 
the evaluation of the noise associated with the above ob- 
servable. As mentioned above, both instrumental noise 
and cosmic variance are considered. We then move to 



8 



estimate the variance of our correlator 

a 2 = (ST 2r K 2 ) - (5F r K ) 2 . (22) 

Since {5T r n) 2 is just the square of the signal, we aim 
here to obtain estimates for {8F 2r K 2 ). From Eq. we 
get: 



XF , ^ l (xc,Xf) 
a{Xc) 



; XF , , w l ( X 'c,xf) [ xq , 



XQ 



'Xi J Xi 

x (<r(^x,)« r (^x;)«(^xc)^,x / c)) > (23) 

where there are now two integrals running along the con- 
vergence los (on Xc and x'c) and two running along the 
Lyman-Q spectrum (on Xq and Xq)- The correlator ap- 
pearing in the integrand of Eq. (|23|) is characterized by 
an even (2r + 2) number of S factors. This implies that an 
approximation to its value can be obtained using Wick's 
theorem. When Wick's theorem is applied, many differ- 
ent terms will in general appear. Adopting for sake of 
brevity the notation S(n,Xi) = terms characterized 
by the contraction of 5i and Sj will receive non-negligible 
contributions over the overlap of the respective los. The 
terms providing the largest contribution to (SF n ) are 
the ones where S c is contracted with 6 C '' these terms 
in fact contain the value of the cosmic variance of the 
convergence and receive significant contributions from all 
points along the los from the observer all the way to the 
last scattering surface. On the other hand, whenever we 
consider the cross-correlation between a 5 C and a 5 q , this 
will acquire a non-negligible value only for those set of 
points where the los to the last scattering surface over- 
laps with the Lyman-a spectrum. As such, these terms 
are only proportional to the length of the Lyman-a spec- 
trum, and thus sensibly smaller than the ones containing 
the variance of the convergence. We note in passing that 
the same argument should also apply to the connected 
part of the correlator, which should be significantly non- 
zero only along the Lyman-a spectrum. Mathematically, 
these facts become apparent from Eq. (f!?3")) above, where 
terms containing (6 c 5 c i) are the only ones for which the 
integration over Xc and x'c can be traded for an integra- 
tion over Axe and an integration over Xc that extends all 
the way to xf- If on the other hand S c is contracted with 
a S q factor, then the approximation scheme of Eqs. (JSJH]) 
leads to an integral over Ax and to an integral over Xq 
that extends only over the length probed by the Lyman- 
a spectrum. It seems therefore possible to safely neglect 
terms where the 5's referring to the convergence are not 
contracted with each other. 



The variance of 5 J- k 

We start by considering the variance of 6 T n. Setting 
r = 1 in Eq. and using Wick's theorem we obtain 

(S q 6 q ,6 c 5 c >) « 2{6 q 6 c ){6 q ,5 c ,) + {6 q 6 q >){6 c 5 c >). (24) 



We notice immediately that the first term is twice the 
square of (6 J 7 k) , while the second term is proportional to 
two correlation function characterized by cutoffs acting 
either on the modes that are parallel or perpendicular to 
the los, but not on both. It is then possible to show that 



(Sqd q > 



D(x q )D(x q >) H^(A Xq ;oo,k L /V2 



4 0) (Ax g ;oo, k/y/2 



(25) 



(8 c 8 d ) = D( Xc )D(xc')H^(Ax c - 7 k c /V2 7 ^) 7 (26) 
(S q 5 c ) = D( Xc )D(xq) 

x [H°(Ax;kc,k L )-H°(A X ;kcM , (27) 



(S 2 ) = D 2 ( X q) H^(x q ,k L )+H^( Xq ,k) ,(28) 



(0), 



where the last two equations have been added here for 
sake of completeness, as they will be useful in what fol- 
lows. The variance of SJ-'k is then 



(5T k) 



A(3 



3H 2 fl T 



2c 2 

x I d Xq D 2 (x q ) 



XF 



dx c ^^D 2 (x c ) 



a 2 (Xc) 

r AXc -° fn\ 



/ dA Xq 

'- A X,,o 

r(o) 



dAxcK>(Ax c ;k c /V2,^) 
H^'(A Xq ;oc, k L /V2) 



-Axc.o 
(0), 



^ 0) (Ax ? ;oo, k/V2) 



(29) 



In the upper panels of Fig. [3] we show the values ob- 
tained for the standard deviation of 8J-K, for two different 
CMB experiments' resolution, again turning off the IGM 
physics evolution and focusing on the growth of struc- 
ture. 



The variance of 8T 2 k 



Setting r = 2 in Eq. (|23j) . we then apply Wick's the- 
orem to (S 2 5 2 ,S c 8 c i). Neglecting again terms where the 
<5 c 's are not contracted with one another, we obtain 



(5 2 5 2 8 C S C ,) « 2{S 2 q 5 c ){S 2 6 c ,) 

+ (S C 6 C ,) ((S 2 q )(5 2 ) + 2(5 q 5 q ,) 2 ), (30) 
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FIG. 3: Estimates of the standard deviation of the correlator &Tk, (upper panels) and 8T 2 k (lower panels) along a single line-of- 
sight as a function of the source redshift z and of the length of the measured spectrum Az, for convergence maps recontructed 
from Planck (left panels) and ACTPOL (right panels). As before, we set A = j3 = 1, effectively turning off the physics of IGM, 
to make apparent the physics of structure formation. 



which then leads to the expression for a\ 



(5T 2 k) 2 



A/3 



2c 2 



Yf 



W2(Xc,Xr) 



a 2 (Xc) 

/ A Xc,0 
dA Xc H { °\A Xc ;k c /V2,oo) 



4°\ Xq ,k L ) + H^( Xq ,k) 

Ax<j,o 



r(o), 



XQ 



dx q D 2 { Xq ) 



9 / n4 



dx q D\ Xq ) I dA Xq H^(A Xq ;^,k L 

i^ 0) (A X? ;oo, k/V2) 2 



(31) 



In the lower panels of Fig. [3] we show the estimates 
for the standard deviation SJ- 2 k along a single line-of- 
sight for the two different CMB experiment. We note in 
Fig. [3] the same trends that have been pointed out for the 



correlator itself in Fig.[T]and[21 the standard deviation of 
bTn and of b~T 2 k increase almost linearly with increasing 
length of the Lyman-a spectrum Az and it decreases as 
the source redshift z is increased because of the fact that 
the spectrum probes regions that are less clumpy. Also, 
by increasing the resolution of the CMB experiment used 
to reconstruct the convergence map, the deviation of 5Thl 
and SJ- 2 k also increase: if on one hand more modes carry 
more information, on the other hand they also carry more 
cosmic variance. 



One last aspect to note here is that while the signal 
for (8T 2 k) arises from a three point correlation func- 
tion (which in the gaussian approximation would yield 
zero), the dominant terms contributing to its variance 
arise from products of two point correlation functions. In 
particular, it is possible to show that the terms appear- 
ing in the second line of Eq. (|30[) significantly outweight 
the square of the signal that appears in the first line. 



10 



0.7 
0.6 

0.5 
0.4 
0.3 
0.2 L 
0.1 h 



-0.050- 



-0.025 - 



N 
<l 



0.7 
0.6 

0.5 

0.4 

0.3 







0.^0-^1 






0.12 5 -"1 






.0.100 -i 






n 075 i 






.0 050 : 


: . . _^ o.n?^_: 



1.5 



2.0 



2.5 



3.0 



1.5 



2.0 



2.5 



3.0 




FIG. 4: Estimates for the signal-to-noise ratios for the observation of the correlators (5J-k) along a single line-of-sight as a 
function of the source redshift z and of the length of the measured spectrum Az, for Planck (left panels) and ACTPOL (right 
panels). As long as the functions A and /3 can be assumed to be constant in the redshift range spanned by the Lyman-a 
spectrum, these result do not depend on the specific value taken by the latter. 



D. Signal-to-Noise ratio 

We now have all the pieces to assess to what extent 
the (SJ- r K) correlations will be detectable by future ob- 
servational programs. Even before moving to plot the 
S/N ratios for 8Tn and 8T 2 n it is possible to point out 
a couple of features of these ratios. First, we note that 
the S/N ratio for bTn and 8T 2 k do present a radical dif- 
ference in their dependence on the QSO source redshift. 
This is because the signal for 5T 2 n is characterized by 
mode coupling, whereas the dominant contributions to 
the variance are not. Physically, the signal for b~T 2 n is 
more sensitive to the growth of structure with respect to 
its variance: while for the former the growth of long wave- 
length modes enhances the growth of structure on small 
scales, for the latter long and short wavelength modes 
grow independently at the same rate. Mathematically, 
this is apparent when comparing Eq. ^ with Eq. pip : 
while the (SJ- 2 k) signal carries four powers of the growth 
factor, the dominant terms contributing to its variance 



carry only six. In this case then the S /N is characterized 
by four growth factors in the numerator and only three in 
the denominator, thus leading to a "linear" dependence 
of S/N on the redshift (modulo integration over the los 
and behaviour of the lensing window function). Note 
that this is in stark contrast with the (SJ-'k) case, where 
the signal is not characterized by mode coupling and the 
number of growth factors are equal for the signal and its 
standard deviation, thus leading to a S/N ratio with no 
dependence on the source's redshift. 

Second, we note that S/N does not depend on the value 
of any constant. In particular, regardless of their redshift 
dependence, the S/N ratio will not depend on the func- 
tions A and (3 used to describe the IGM. This is of course 
very important since in such a way, at least in linear the- 
ory and using the FGPA at first order, the dependence 
on the physics of the IGM cancels out when computing 
the S/N ratio. 

In Fig. [3] we show the estimates for the S/N per los of 
the (SJ-k) (upper panels) and (SJ- 2 k) (lower panels) mea- 
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CMB Exp. 


S/N 


Total S/N 


Total S/N 


CMB Exp. 


S/N Total S/N 


Total S/N 




per los 


in BOSS 


in BigBOSS 




per los in BOSS 


in BigBOSS 


Planck 


0.075 


30 


75 


Planck 


0.024 9.6 


24 


ACTPOL 


0.130 


52 


130 


ACTPOL 


0.05 20.0 


50 



TABLE I: Estimates of the total and per single los signal-to- 
noise (S/N) of the {5 J- k) cross-correlation for different CMB 
experiments combined with BOSS and BigBOSS. 



TABLE II: Estimates of the total and per single los signal- 
to-noise (S/N) of the {5T 2 k) cross-correlation for different 
CMB experiments combined with BOSS and BigBOSS. 



surements. As expected, while the S/N for (8!Fk) does 
not show any strong redshift dependence, the S/N for 
{8T 2 k) decreases linearly with increasing source redshift: 
the growth of structure is indeed playing a role and shows 
that QSOs lying at lower redshift will yield a larger S/N. 
Also, in both cases an increase in the resolution of the 
experiment measuring the convergence field translates in 
a larger S/N and in a larger derivative of the S/N with 
respect to Az. This is not surprising, as it is reasonable 
to expect that a higher resolution convergence map will 
be carrying a larger amount of information about the 
density field. 

All this suggests that depending on what is the corre- 
lator that one is interested in measuring, different strate- 
gies should be pursued. In case of (SFk) increasing the 
length of the spectra will provide a better S/N. In case 
of (8T 2 k), however, Fig. (01 suggests that an increase 
in the number of quasar will be more effective in pro- 
ducing a large S/N, whereas an increase in the redshift 
range spanned by the spectrum will increase the S /N only 
marginally. 

Having obtained the S/N per los, we can then estimate 
the total S/N that will be obtained by cross-correlating 
the BOSS sample (1.6 • 10 5 QSOs) and the proposed Big- 
BOSS sample [l0[ (10 6 QSOs) with the convergence map 
measured by Planck or by the proposed ACTPOL ex- 
periment considered. Assuming a mean QSO redshift 
of z = 2.5 and a mean Lyman-a spectrum length of 
Az = 0.5, a rough estimate of the S/N for the measure- 
ments of {8Tn) and of {8T 2 n) are given in Tab.HlandHIl 



It is necessary to point out here that despite that 
the value of the S/N for (8Tk) is almost three times 
larger than the one for (8J- 2 n), the actual measurement 
of the former correlator strongly depends on the ability 
of fitting the continuum of the Lyman-a spectrum. The 
(<5.F 2 re) correlator, on the other hand, is sensitive to the 
interplay between long and short wavelength modes and 
as such should be less sensitive to the continuum fitting 
procedure. Therefore, even if it is characterized by a 
lower S/N, it may actually be the easier to measure in 
practice. The numbers obtained above are particularly 
encouraging since the S/N values are typically very large 
and well above unity. 



E. Analysis 

Having developed a calculation framework for estimat- 
ing (8!F r K) and the S/N for their measurement, we turn 
to estimate what is the range of Lyman-a wavelengths 
contributing to the signal and what is the effect of chang- 
ing the parameters that control the experiments' resolu- 
tion. 



Spectral Analysis 

We investigate here how the different Lyman-a modes 
contribute to the correlators. This should tell us whether 
long wavelength modes have any appreciable effect on 
our observables and what is the impact of short and very 
short wavelength modes (in particular the ones that are 
expected to have entered the non- linear regime). 

Since the mean flux T appearing in the definition of 
the flux fluctuation 8 J- = ( T — T) jT is a global quantity 
which is usually estimated from a statistically significant 
sample of high resolution QSO spectra (see the discus- 



sion in [44J for the impact that such quantity has on 
some derived cosmological parameters), 8T is sensitive 
also to modes with wavelengths longer than the Lyman- 
a spectrum. These modes appear as a "background" in 
each spectra but they still have to be accounted for when 
crosscorrelating 8 J- with re because the fluctuation in the 
flux is affected by them. More specifically, a QSO that 
is sitting in an overdense region that extends beyond the 
redshift range spanned by its spectrum will see its flux 
decremented by a factor that in its spectrum will appear 
as constant decrement. On the other hand, if the QSO 
spectrum extends beyond the edge of such overdensity, 
this mode would appear as a fluctuation (and not as a 
background) in the spectrum. This extreme scenario is 
somewhat mitigated by the fact that present and future 
QSO surveys will have many QSOs with los separated by 
few comoving Mpc [4fj| : as such, fluxes from neighbor- 
ing QSO lying in large overdense regions should present 
similarities that should in principle allow to detect such 
large overdensities in 3D tomographical studies [30l ]. 

To measure the contributions of the different modes to 
the correlators, we vary k\ and fcz, to build appropriate 
filters. As can be seen from Fig. where three such 
filters are plotted for {fc; = 0.001, k L = 0.01}, {h = 
0.01, kr, = 0.1} and {ki = 0.1, kr, = 1}, the gaussian 
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i7 



k'L 



1(^)1 



Ratio 



1.00e-04 1.00e-03 1.66e-04 1.77e-04 9.39e-01 
1.00e-03 1.00e-02 1.20e-03 1.21e-03 9.87e-01 
1.00e-02 1.00e-01 2.12e-04 6.29e-04 3.37e-01 
1.00e-01 1.00e+00 6.11e-07 1.42e-03 4.30e-04 
1.00e+00 1.00e+01 7.26e-08 2.44e-03 2.97e-06 

TABLE III: Contribution of the different wavenumbers (split 
over decades) to the absolute value of the correlator (STu), 
its standard deviation a&Ftt and ratio of the two quantities. 
In this calculation we took into account the evolution of A 
with redshift. 



A". 



Ratio 



1.00e-04 1.00e-03 1.08e-04 2.18e-02 4.99e-03 

1.00e-03 1.00e-02 6.69e-03 1.96e-01 3.40e-02 

1.00e-02 1.00e-01 5.92e-02 1.31e+00 4.52e-02 

1.00e-01 1.00e+00 3.39e-01 7.06e+00 4.80e-02 

1.00e+00 1.00e+01 9.92e-01 2.07e+01 4.79e-02 

TABLE IV: Contribution of the different wavenumbers (split 
over decades) to the correlator {S^k}, its standard deviation 
<Jst 2 k and ratio of the two quantities. In this calculation we 
took into account the evolution of A with redshift. 



functional form assumed for the window function does 
not provide very sharp filters (hence this spectral analysis 
will not reach high resolution). Also, if = 10 ki then 
the filters add exactly to one. This allows us to measure 
the contributions of the different wavenumber decades to 
the correlators and its standard deviation. 

Tabic Mil and IIVI summarize the results for (8!Fk) and 
(8J- 2 k) respectively. Considering (SJ-k) wc note immedi- 
ately that the signal and the S/N ratio both peaks around 
k ~ 10" 2 /iMpc -1 , as expected from the fact that this 
signal is proportional to the two point correlation func- 
tion, which in turn receives its largest contribution from 
the wavelengths that dominate the power spectrum: iso- 
lating the long wavelength modes of the Lyman-a flux 
would allow to increase the S/N. However, this procedure 
is sensibly complicated by the continuum fitting proce- 
dures that are needed to correctly reproduce the long 



1.0 




1.0E-O4 1.0E-03 1.0E-02 1.0E-01 1.0E+00 1.0E+01 

k (h Mpc"') 

FIG. 5: Three filters used to calculate the contribution of the 
different modes to the correlators, their variance and the SN 
ratio. The filters have {ki — 10~ 3 ,fc.L = 10~ 2 } (solid curve), 
{ki = 10~ 2 ,k L = KT 1 } (dotted curve) and {k t = 10 -1 ,ftz, = 
1} (dashed curve). Also shown is the sum of the filters (red 
dashed-dotted curve). 



wavelength fluctuations of the Lyman-a flux. The be- 
havior of the variance is interesting, as in the first three 
decades shows an oscillating behavior. This is due to the 
different weights of the two terms appearing in Eq. ([29]) 
for each range of wavelengths. In particular, for k < 10~ 2 
/iMpc -1 the variance of (SFn) is dominated by the first 
term, that is just the square of the signal. However, as 
the signal gets smaller with increasing k, for k > 10" 1 
/iMpc" 1 it is the second term that dominates the vari- 
ance. 

Regarding (5T 2 k), it is necessary to point out two 
aspects. First, short wavelengths (high-fc) modes pro- 
vide the larger contribution to both the correlator and 
its standard deviation. Second, for k > 10~ 2 /iMpc -1 
the ratio of the contribution to the correlator and to its 
standard deviation remain almost constant. This means 
that above 10~ 2 h Mpc" 1 the different frequency ranges 
contribute roughly in the same proportion. This fact is 
both good news and bad news at the same time. It is 
bad news because it means that increasing the resolution 
of the Lyman-a spectra does not automatically translate 
into increasing the precision with which the correlator 
will be measured, as the high-/c modes that are intro- 
duced will boost both the correlator and its variance in 
the same way. On the other, this appears also to be 
good news because it tells us that low resolution spectra 
which do not record non-linearities on small scales can 
be successfully used to measure this correlation. To in- 
crease the S/N ratio and to achieve a better precision for 
this measurement it is better to increase the number of 
QSO spectra than to increase the resolution of each single 
spectra. Finally, cutting off the long-wavelength modes 
with k < 10" 2 h Mpc" 1 should not have a great impact 
on the S/N ratio or on the measured value of the corre- 
lator: if on one hand the contribution of the modes with 
k < 10~ 2 /iMpc -1 are noisier due to cosmic variance, on 
the other hand the absolute value of such contributions to 
the correlator and to its variance are negligible compared 
to the ones arising from k > 10~ 2 /iMpc" 1 . We can see 
this fact also comparing the last column of Tab. IIVI with 
the right panel of Fig. [S] where the absolute value of the 
S/N ratio is plotted for varying values of the cutoffs fcj, 
and kc- By looking at the last column of Tab. IIVI we see 
that the ratio between the correlator and its standard de- 
viation increases until about k ~ 10 _2 /iMpc" 1 where it 
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FIG. 6: Value of {ST 2 ^) (left panel, red dashed contours), of its standard deviation (left panel, black solid contour) and of 
its S/N ratio (right) for a single QSO lying at z = 2.6 and whose spectrum covers Az — 0.5. Here we assume ki = and 
A = p=l. 



levels off. Looking at the right panel of Fig. [5] we notice 
exactly the same trend: increasing the resolution of the 
spectrum kt above 10 _2 /iMpc _1 does not improve the 
dramatically the S/N ratio. This is because from that 
point on each new mode contributes in almost the same 
amount to the correlator and to its standard deviation. 



Dependence on experimental resolutions 



To analyze the impact of a change in the resolution of 
the experiments measuring the CMB convergence map or 
the Lyman-cv flux we consider a single QSO at rcdshift 
zq = 2.6 whose spectrum covers Az = 0.5 and vary kL 
and kc- In this case we set ki = 0. 

In Fig. [5] we show the value of ((5JF 2 k), of its standard 
deviation and of its S /N ratio for varying values of fc^ and 
kc- We note that both the correlator and its standard 
deviation increase with increasing resolution: this makes 
physical sense as increasing the resolution increases both 
the amount of information carried by each experiment 
and the cosmic variance associated with it. Except for 
very low values of kc, an increase in the resolution of the 
Lyman-cv spectrum is characterized by an almost equal 
amount of increase in both the correlator and its cosmic 
variance. This implies that the S/N becomes roughly 
constant for kL ^ 10~ 2 ft,Mpc~ . On the other hand, 
increasing kc increases both the correlator and its cosmic 
variance only up to the point where kc — kL- 



III. COSMOLOGICAL APPLICATIONS 
A. Neutrinos 

Massive neutrinos are known to suppress the growth of 
structure in the early universe on intermediate to small 
scales k > 1CT 2 /iMpc" 1 [II]. Since (&F 2 /c) is mostly 
sensitive to the same range of scales, it seems reasonable 
to examine to what extent massive neutrinos will alter 
the (S^k) signal. The argument could also be turned 
around, asking how well a measurement of (8J- 2 k) would 
allow to constrain the sum of the neutrino masses. In 
this first work, we take the first route and we simply 
calculate how the (ST 2 ^) signal is affected by different 
values of the neutrino masses. We leave the analysis of 
the constraining power of (S^n) to a forthcoming work. 

Quite generally massive neutrinos affect the matter 
density power spectrum in a scale dependent way (see 
p6| for a review). To account for this effect in an ex- 
act way it would require substantial modifications of the 
formalism and of the code that we are currently using 
to evaluate (5T 2 k). In particular, it would not be pos- 
sible any longer to separate the integrations over the co- 
moving distance from the ones over the wavenumbers k. 
We leave this important development to a future project 
and for the purpose of this work we rely on the follow- 
ing approximation (47j for the growth of the dark matter 
perturbations 



6 cdm cx Dia) 1 i f " 



(32) 



where /„ = f2„/Q m . The above expression should be ac- 
curate from the very large scales down to those mildly 
non- linear ones of the Lyman-cv forest. Departures at 
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FIG. 7: Cross-correlation of the variance of the Lyman-a flux and CMB convergence as a function of redshift for the three 
different cosmological models with massive neutrinos shown in Table IV. Black, orange and red (with la) error bars refer 
to Em„(el/) = 0.54,0.4,0.15, respectively. The black dashed line shows the prediction for a massless neutrino cosmology 
consistent with WMAP-5 data. Four different cases are reported here for Planck+BOSS (top left), Planck + BigBOSS (bottom 
left), ACTPOL+BOSS (top right) and ACTPOL+BigBOSS (bottom right). The redshift evolution of A is here taken into 
account. 



small scales are best handled with N-body or hydrody- 
namical codes [H|. 

The second aspect that we need to take into account 
before proceeding with the calculation is that consistency 
with CMB data requires that a change in the sum of the 
neutrino masses is accompanied by a change in the power 
spectrum normalization as [HJ . This fact has a profound 
consequence. Just by counting the number of powers of 
the power spectrum that enter in the different expres- 
sions, it is straightforward to note that (S^k) ~ erf, 
that c^5jr2,,) ~ tTg and that its S/N ratio is proportional 
to erg- Consequently, a change in the neutrino masses, 
which requires a change in erg to maintain consistency 
with CMB data, will cause a change in (5J- 2 k). 

To take this into account we proceed as follows. First 
we consider the set of values allowed by the WMAP-5 
data in the og — Sm„ space at 95% CL. These corre- 
spond the the dark red area of the center panel of Fig. 
17 in Komatsu et al. 1121. We then choose three flat mod- 



Num. 






Q v 


£m„(eV) 


178 


h 


1 


0.269 


0.719 


1.2e-2 


0.54 


0.657 


0.70 


2 


0.269 


0.722 


8.8e-3 


0.40 


0.708 


0.70 


3 


0.269 


0.728 


3.3e-3 


0.15 


0.786 


0.70 


4 


0.256 


0.744 


0.0 


0.0 


0.841 


0.72 



TABLE V: Values of the cosmological parameters assumed 
to estimate the effect of massive neutrinos on {ST «). All 
models assume flat geometry. 



els with massive neutrinos consistent with the WMAP-5 
data and we use CAMB to generate the respective dark 
matter power spectra to be used in the calculation. The 
value of the cosmological parameters used for each model 
are summarized in Tab. [V] 

One last point is left to be considered. Note in fact 
that the S/N ratio for (SJ^k), although increasing with 
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0.1 



1.0 



FIG. 8: Growth factors for the WMAP-5 fiat ACDM cos- 
mology (dotted curve) and for the early dark energy (EDE) 
model assumed in this section for comparison (solid curve). 



Az, does not increase at a very high rate. It seems there- 
fore possible to speculate that subdividing the Lyman-a 
spectra into sub-spectra, each of length dz — 0.1, de- 
spite lowering the S/N ratio for each single sub-spectrum, 
would allow to reach a better measurement of the rcdshift 
dependence of the signal. 

Figure [7] below shows the result of applying the lat- 
ter procedure. The black, orange and red data points 
represents predicted values of the (8!F 2 k) correlator for 
values of m v = {0.54,0.4,0.15} respectively, while 
the dashed black line shows the value of the correlator 
for a ACDM cosmology with massless neutrinos. As one 
can see, the cross-correlation signal is quite sensitive to 
the presence of massive neutrinos and already BOSS and 
Planck could provide constraints on the strength of such 
correlators. As pointed out above, this is due to the fact 
that more massive neutrinos requires smaller values of 
us, which in turn depresses the signal. 

It is here necessary to point out one important caveat. 
In this paper, we are making a tree-level approximation 
to the growth rate of k modes: this enables us to separate 
integrations along the comoving distances from integra- 
tions on the different modes. As previously mentioned, 
this approximation does not include the scale-dependent 
effects of neutrinos on the growth rate of structure. Sim- 
ilarly, this approximation also does not allow us to take 
into account the non-linearities induced by gravitational 
collapse, which on the other hand tend to enhance the 
power spectrum on small scales. We will need to either 
use Hyper-Extended Perturbation Theory results or non- 
linear simulations to evaluate these effects. 



QSO sample 


CMB Experiment 


Ax 2 


1.6 ■ 10 5 (BOSS) 


Planck 


0.3451 


1.6 ■ 10 s (BOSS) 


ACTPOL 


2.157 


1.0 • 10 6 (BigBOSS) 


Planck 


1.458 


1.0 ■ 10 6 (BigBOSS) 


ACTPOL 


9.117 



TABLE VI: Summary of the estimated Ax 2 between EDE 
and ACDM for four different combinations of future QSO and 
CMB experiments using the (5J- 2 k) correlator. 



B. Early Dark Energy 

Since early dark energy or deviations from general rel- 
ativity affect the growth rate of structure as a function 
of scale, the measurements of (SJ 72 k)(z) can in principle 
probe these effects. Here we focus on early dark energy 
(EDE) models, where dark energy makes a significant 
contribution to the energy density of the universe over 
a wide range of redshifts. The differences between EDE 
models and pure ACDM are particularly evident at high 
redshifts, when the former has been shown to influence 
the growth of the first cosmic structures both in the linear 
and in the non-linear regime. 

We consider here the EDE model proposed in [pi and 
recently constrained by 50] (model EDE1 of [50]). We 
compare this model with the ACDM cosmology assumed 
until now. The differences in the growth factors for these 
two models is shown in Fig. [5] (the difference in the Hub- 
ble parameter evolution is smaller). 

We quantify the departure of the correlators predicted 
for the EDE model from the ACDM one using the fol- 
lowing expression: 



AX 2 = E 



((^k)ede - (&F"«)acdm) 

2 



(33) 



The results are shown in Fig. [5] and are summarized 
in Table IV. In this case the differences between EDE 
and ACDM are very limited and could only be appre- 
ciated at some significance with an advanced CMB ex- 
periment like ACTPOL and by increasing the number of 
spectroscopic QSOs with BigBOSS. However, it is worth 
stressing that the two models presented here are in per- 
fect agreement with all the low redshift probes and the 
large-scale structure measurements provided by galaxy 
power spectra, CMB, Type la supernovae and Lyman-a 
forest. Therefore, possible departures from ACDM can 
be investigated only exploiting the capabilities of this in- 
termediate rcdshift regime with such correlations or with 
similar observables in this redshift range. 



IV. CONCLUSIONS 

This work presents a detailed investigation of the cross- 
correlation signals between transmitted Lyman-a flux 
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FIG. 9: Value of (S^k) estimated for the early dark energy EDE model of [4jj. The dashed black line shows the expected 
value of the correlator for the ACDM cosmology assumed so far. Four different cases are reported here for Planck+BOSS (top 
left), Planck + BigBOSS (bottom left), ACTPOL+BOSS (top right) and ACTPOL+BigBOSS (bottom right). The redshift 
evolution of A is here taken into account. 



and the weak lensing convergence of the CMB along 
the same line-of-sight. One of the motivations behind 
this work is that the Lyman-a forest has already been 
shown to be a powerful cosmological tool and novel 
ways of exploring and deepening the understanding of 
the flux/matter relation could significantly improve our 
knowledge of the high redshift universe. These corre- 
lators are able to provide astrophysical and cosmolog- 
ical information: since they are sensitive to both the 
flux/matter relation and the value of cosmological pa- 
rameters, in principle they can be used to put constraints 
on both. 

The correlators investigated in the present work have 
a clear physical meaning. The correlation of 8T with k 
measures to what extent the fluctuations along the los 
mapped by the Lyman-a forest contribute to the CMB 
convergence field. This correlation is dominated by long 
wavelength modes (fc < lO^ 1 /iMpc -1 ) and as such is 
more sensitive to Lyman-a forest continuum fitting pro- 
cedures. The correlation of the flux variance 8T 2 with k 



measures to what extent the growth of short wavelength 
modes (mapped by the Lyman-a flux) is enhanced or de- 
pressed by the fact that the latter are sitting in regions 
that are overdense or underdense on large scales. This 
interplay between short and long wavelength modes is 
well exemplified by the redshift dependence of the S/N 
ratio for (8!F 2 k): lowering the redshift increases the S/N 
ratio because while the variance of (8T 2 k) is dominated 
by the independent growth of long and short wavelength 
modes, the value of (8T 2 n) itself receives an extra con- 
tribution due to the fact that the growth of the short 
wavelength modes is enhanced by the presence (and in- 
dependent growth) of the long wavelength modes. Fur- 
thermore, this correlator is sensitive to intermediate-to- 
small scales (k > 10 -2 /iMpc -1 ) and as such it should be 
less sensitive to Lyman-a forest continuum fitting proce- 
dures. 

To estimate the values of the correlators, their variance 
and their S/N ratio we rely on linear theory and simple 
approximations, such as the fluctuating Gunn-Pcterson 
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approximation at first order. Although the framework is 
simplified, the results are by no means obvious since dif- 
ferent modes enter non-trivially in these quantities and 
in their signal-to-noise ratio. We estimate that such cor- 
relations may be detectable at a high significance level by 
Planck and the SDSS-111 BOSS survey, experiments that 
are already collecting data. Moreover, our investigation 
of the modes of the Lyman-a forest that contribute to 
(S^k) shows that the low-resolution Lyman-a spectra 
measured by SDSS-111 (which is aimed at the measure- 
ment of BAO at z = 2 - 4 [13, [5l|) should have enough 
resolution to yield a significant S/N. 

The peculiar dependence of {8!F 2 k) on intermediate- 
to-short scales and its sensitivity to the value of the power 
spectrum normalization a% makes it a very useful cosmo- 
logical tool to test all models characterized by variations 
of the power spectrum on such scales. In particular, 
we applied our estimates to evaluate the sensitivity of 
(8T 2 k) to changes in og due to variations in the sum 
of the neutrino masses and to show how promising this 
measurement could be in constraining the latter. 

Finally, some caveats are in order. First, the code de- 
veloped to estimate (5J- 2 k) and its variance is based on 
the tree-level perturbation theory results reported here. 
As such, the results shown do not take into account non- 
linearities induced by gravitational collapse. The ex- 
tension of the analytic results to take into account this 
aspect is actually quite straightforward, as it only re- 
quires the implementation of the so-called "HyperEx- 
tended Perturbation Theory" for the bispectrum [4p| . 
However, the implementation of such changes in a nu- 
merical code are less trivial, as the integrations over the 
power spectrum and over the comoving distance cannot 
be factored any longer. We have nonetheless reason to 
speculate that the nonlinearities induced by gravitational 
collapse will not dramatically change the picture outlined 
here. At the redshift range spanned by the Lyman-a for- 
est nonlinearities are normally mild and confined to short 
scales. Furthermore, as shown in Sec. Ill El the S/N ratio 
for {5T 2 k) dominated by modes with k > 10~ 2 ftMpc -1 , 
but all decades above 10~ 2 /iMpc -1 contribute in the 
same proportion to both the signal and its variance. It is 
therefore conceivable to filter out of the Lyman-a spec- 
tra the shortest scales, which arc the most affected by 
nonlinearities and still be able to retain a non-negligible 
S/N. 

The second caveat pertains the estimate of the correla- 
tors' variance. It is in fact necessary to point out that to 
obtain such estimates Wick's theorem has been applied. 
Whether the use of Wick's theorem may or may not lead 
to an accurate result when considering the variance of 
(5J- 2 k) is debatable. On one hand it is possible to point 
out that the largest part of the signal arises at small sep- 
arations, where the value of the correlator is dominated 
by its connected part. Analogously, it could be possi- 



ble to argue that the use of Wick's theorem may lead to 
underestimating the correlators' variance. An exact eval- 
uation of the variance of {5!F 2 k), however, requires the 
exact calculation of a six point function, that to the best 
of our knowledge has never been determined. On the 
other hand it is also possible to point out that the con- 
nected part of (S C 5 C 'S 2 S 2 ,) will be significantly non-zero 
only when the distances between the different points is 
small. As such, this term will give a non-zero contribu- 
tion proportional to the length of the Lyman-a spectrum, 
which should be subdominant with respect to the ones 
considered in section III CI that are proportional to the 
distance from the observer all the way to the last scat- 
tering surface. 

The third caveat pertains the expansion of the expres- 
sion for the flux, Eq. (JTJ) . Despite the fact that the ex- 
pansion carried out in Eq. ^ is correct on scales larger 
than about 1 /i Mpc, we point out here that the flux as 
expressed in Eq. ([T]) is intrinsically a non-linear function 
of the overdensity field. It is therefore reasonable to won- 
der whether the non-linearities induced by this non-linear 
mapping would somehow affect the conclusions presented 
here. A simple way to sidestep the present question is 
to undo the non-linear mapping by defining a new ob- 
servable T = — ln(.F) = A(l + (Sigm)' 3 and proceed by 
measuring its correlations. 

The best way to assess to what extent the above 
caveats affect the estimates reported in the present work 
is through numerical simulations, calculating the conver- 
gence field on a light cone and at the same time measuring 
Lyman-a forest synthetic spectra and cross-correlating 
the two. This will be the next step in our investigation 
and the focus of the next publication. 

Finally, on the analytical side we still need to address 
the estimate of the correlators when the power spec- 
trum shows evolution in redshift and on different scales 
at the same time. As pointed out, (5T 2 k) is sensitive 
to scales k > 10 -2 /iMpc -1 . As such this correlator is 
an ideal tool to test modifications of gravity that show 
scale dependent growth. At the same time, this devel- 
opment would also allow the implementation of the hy- 
perextended perturbation theory results and as such to 
address analytically the impact of gravity induced non- 
linearities on the value of the correlators. 
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APPENDIX A: DERIVATION OF PERTURB ATIVE RESULTS FOR {8T 2 k) 

In this section we derive the expression for (8T 2 k) shown in the text, Eqs. (|19H21j) . We move from Eq. ([5]) and need 
to find an efficient way to evaluate (S 2 (n, x g )i5(n, Xc))- We start by Fourier transforming this cumulant correlator to 
get 

(S 2 q S c ) a (k2,\\)W K (h tl _)(S(h)8(k2)S(h)) 

^ + fc3)^ Q (fc 1J |)^ Q (fc 2i ||)^ K (fc3,x)B(fcl,fc2,fc 3 ) 



F 2 (k 1 ,k 2 )P L (k 1 , X i) P L (k 2 ,X2) + F 2 (k 2 ,k 3 )P L {k 2 ,X2) P L (k 3 ,X3) + F2(k 3 ,k 1 )P L (k 3 ,x 3 ) PL(ki, X i) ■ 

(Al) 



d 3 k x 


d 3 k 2 


d% 


(2tt) 3 


(2n) 3 


(2ir) 3 


d 3 h 


d 3 k 2 


d% 


(2tt) 3 


(2tt) 3 


(2tt) 3 


d 3 h 


d 3 k 2 


d 3 k 3 


(2tt) 3 


(2tt) 3 


(2tt) 3 



x 2 



In the second line we introduced the bispectrum B{k\,k 2l k 3 ), while in the third line we replaced the bispectrum 
with the expression for its kernel F 2 and products of the linear matter power spectrum Pi(fc, For sake of brevity 
we keep implicit the dependence of the window functions on the cutoff scales: W a (k i ^) = W a (k i ^,ki J ,ki) and 

W K (k^x) = W K (k^x,kc)- Next, we point out that the evaluation of Eq. (|A1[) requires in general the integration over 
a six dimensional fe-space, which is further complicated by the fact that the different window functions break the 
spherical symmetry that one would normally exploit. 

In what follows we adopt the tree level approximation to the bispectrum kernel, 

= Uf^ {kf + k > )+ l(^) ' (A2) 
which can readily be obtained from the more general expression derived by Scoccimarro and Couchman (4p| 

^f £PT (^,fc J ) = ^a(n,h)a(n,k J ) + ^^^(kf + kf)b(n,k t )b(n,k 1 ) + ^^^\ c(n, h) c(n, (A3) 



setting the three auxiliary functions a(fc), b(k) and c(k) that allow to account for non- linear growth of structure equal 
to unity. A generalization of the results shown below to take into account the more general formulation of Eq. (|A3[) 
is straightforward to derive. 

To proceed further we note that each of the three terms appearing in the square bracket of Eq. (|A1[) depend only 
on two of the three wavevectors. When moving from the second to the third line, it is then essential not to carry out 
the integration over the delta function, because for each of these terms we integrate the Dirac S in order to obtain an 
expression that depends only on the same wavevectors that appear in the F 2 kernel. The fact that two of the three 
physical points are the same also spoils the cyclic symmetry of the bispectrum. In particular, the {1,2} term will 
differ from the {2, 3} and {3, 1} terms. We therefore let 

(S 2 q S c ) = (S 2 q S c ) 1 . 2 + 2(S 2 5 c ) 2 , 3 , (A4) 

and start by considering (S^S^i^. Integrating over the Sd function in order to get rid of &3 in favor of k\ and k 2 , and 
then adopting a cylindrical coordinate system in fc-space we get 

d6 ± F 2 {h ,k 2 )W K [\k h± + k 2< i_ |] • (A5) 



As also recognized in [52j, the most challenging part of the calculation consists of the integration over the angular 
variables. This is because the convergence window function depends on + k 2 .±\. The integration over the 
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angular variables in this case does not necessarily lead to an expression that may be numerically efficient to evaluate. 
In particular we aim to keep integrations factored as much as possible. Our first goal then is to integrate 



# d9_ L F 2 (ki,k 2 )W K [\ki,± + k 2 ,±\} = 2tt cxp 
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K l,. 



k 2 c 



exp 



k c 



d9j_F 2 (ki,k 2 ) exp 



,fcl,J_fc2,_LCOS(0_L) 



(A6) 



where 9± is the angle between k\ t x an d &2,x- Now, as far as the integration over the angular variable is concerned, 
the kernel F 2 can be written as 



F 2 (ki,k 2 ) =R + S cos(6>_l) + T cos 2 



(A7) 



where we have decomposed k into its component parallel and perpendicular to the los according to k = k\\n + k± and 
extracted the terms that are proportional to different powers of cos(9±) 
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Integration over the angular variable can then be carried out by remembering that 
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(A9) 
(A10) 

(All) 
(A12) 
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where /„ denotes the modified Bessel function of the first kind and n-th order. The integration over the angular 
variables yields 



dcj) / d9 ± F 2 (k u k 2 ) W c [\k h± + k 2 , ± \] = (2tt) 2 exp 
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The difficulty with this result is that every term depend on the product fci j j_fc2 j j_. As such, we are facing a 2D joint 
integration over the whole [fci k 2t j_] domain. If on one hand this is doable, on the other hand we are more interested 
in obtaining a final result which is a product of integrals instead of the integral of the product. It is possible to move 
around this obstacle recalling that (Abramowitz and Stegun [53[, 9.6.10) 
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We can write the modified Bessel function splitting the dependence on fci.j_ and k 2 ^± as 
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(A14) 

(A15) 
(A16) 
(A17) 



20 



where for sake of brevity we use the following notation for the coefficients 
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00 
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(A20) 



Now, however complicated, this form allows us to factor the different integrals. Let's start by considering the term 
RIq. We have 
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(A21) 



where in going from the first to the second step we expressed k 2 ^_ as a function of k and k\\ using the fact that 
k 2 = k' 2 + k'j_ . This is necessary because the power spectrum is function of k and not of k± . We can then proceed by 
defining the following functions 
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It is important to notice that because of the integration domain all the above functions are even in ku , regardless of 
the value of m. With the help of these functions we then have 
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We have therefore succeeded in obtaining an expression that has the dependence on k±u and k 2 .\\ completely factored. 
The sums over m and the fact that each term is a product of factors that only depend either on ki ii or on fc 2 .|| allows 
an integration term by term and at the same time to bypass the two dimensional joint integration. 

We can then proceed exactly in the same way for the other two terms, S I± and T I 2 with the only difference that 



in order to obtain expressions where only the coefficients of the modified Bessel function of 0-th order 1q appear we 
use the fact that l[ m ^ = (m + 1) I^ n+1 \ We then obtain for the S term the following expression 
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Finally, the T term gives 



{l,Xl) 



(2nf 



d9±T cos 2 (e±) W c [\ki,± + k 2 ,±\ 



7 



^ m(2m- 1) L m (ki t \\,xi)L m (k2 t \\,Xa.)- 



(A26) 



With the introduction of the definitions (|A22HA23|) and with the series expansion for the modified Bessel function we 
have therefore managed to carry out the integration over the perpendicular part of the wavevector. We are then left 
with the integration over ku . First recall that the window functions acting on the Lyman-a flux are 



W a (k\\,k L ,ki) = 



I _ e -(fc||/fci) 2 



-(fell A) 



(A27) 



and that in Eq. (|A5[) they decouple from one another. We can proceed further by defining the following function 



f£>(A X ,X,kc,k L ) ee / ( ^] exp 



-yl + ^ii A x 



/m(fc||,x;fc c ), 



(A28) 



where / = {H, L}. It is straightforward to note that because all the tilde functions are even in ku, depending 
on the value of n the above Fourier transforms are either purely real (if n is even) or purely imaginary (if n is 
odd). Furthermore, if n is even the above functions are real and even, while if n is odd the above functions are 
imaginary and odd. Carrying out the integration on fcy is then straightforward, as it just corresponds the replacement 

kf\f m (k\\,x\kc) -» fc£/m '(Ax,x;kc,k L ) - k n fffl '(Ax,X\kc,k)- Finally, from a computational point of view this 
approach is rather efficient, as the tilded function need to be calculated only once and then used to construct the 
two-index functions. 

With the help of these auxiliary functions we can finally obtain the following expression for the cumulant correlator 



m=0 



2 E \ | \H^(A x , Xq -k c ,k L )-H^{A x , Xq -k c ,k) 



+ \k L H£\A X , Xq ; k c , k L ) -kH^(A X , Xq ; k c , k)] [k L Lg\A X , Xq ; k c , k L ) - kL$(A X , Xq ; k c , k) 



H£\A X , Xq ; k c , k L ) - H^(A X , Xq ; k c , fc)J [L^(A X , Xq ; k c , k L ) - L%(A X , Xq ; k c , k) 

2 r - - 1 2 

- \ k 2 L L$(A X , Xq ;k c , hi) - P L$(A X , Xq ;k c , k) 

4m 

m(2m — 1) 
7 



•.(»>< 



k L L<£>(A X , Xq ;k c ,k L ) - k L { r ^(A X , Xq ;k c ,k) 



k c 



L<£>(A X , Xq ;k c , k L ) - Lt\A X , Xq ;k c , k) 



(A29) 



A cautionary note is in order. As mentioned above, the functions defined through Eq. (|A28j) are purely imaginary 
if the index (n) is odd. However, notice that in Eq. (|A29j) above there are always two such functions that appear 
together (as in the case with H^L^)), thus ensuring that (<5 2 <5) 12 is always real valued. 

Let's now move to calculate (<5 2 <5)2,3- Notice incidentally that this term is exactly equal to (<5 2 <5) 3^1 . We start from 
the now usual expression 

^M^^^„Ax^ ( _ fc2>|| _ fc3i||)Wa(fc2i||) r- Mb 



<<™>2,3 = 2 



2tt 2tt 

/ dep [ de x F 2 (k 2 ,k 3 )W K (k 3 , x ), 



\k 2 



P(k 2 ,X2) 



k 3 dk 3 
(2tt) 2 



P(k 3 , X3 ) 



(A30) 



where, as previously, we have traded the integrations over ki _± for the ones over ki. In this case the integration over 
the angular variables doesn't pose any problem as the window function W K is actually a function of k 3 ^± only and it 
can be safely pulled out of the angular integrals 



dep / de ± F 2 {k 2 ,k 3 ) 



(2^) 2 



(5 + 1) 



(2nf 



k 



2,11 K 3,|| 



1 



1 



k\ fc| 



(2tt) 5 



h 2 h 2 
, K 2,|| ft 3,|| 

k\ k 2 



k l\\ 
k\ 



k 2 

K 3 ,\\ 

k 3 



iA31) 
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It is here necessary to point out that since W K depends only on A^j., the tildcd functions that will appear when the 
integration over k 2 is carried out will contain no hlter function. We characterize this functions by substituting to kc 
the oo symbol, as to all extent the gaussian filter with kc — > oo just yields unity. We then have 

P(h,X3) [d(f> [ d6 x F 2 {k 2 ,k 3 )W K (k 3 , x ) 



p oo 

^ l fc 2,|| | 


k 2 dk 2 ,? > 

(27r)2 P(fc 2 ,x 2 ) 


p OO 
^ 1^3,11 


k 3 dk 3 

| (27T) 2 


/*OC 
^ |fc2,|| | 


kidk 2 -. 


poo 
'l*a.ll 


k 3 dk 3 
(2tt) 2 


x (2tt) 2 


' 5 1 , , / 

7 + 2 fe2 .ll fe3 .ll 1 


' 1 

V^2 


f 4) 



-P(fe,X3)exp - 



1 1 



P - P 



7 fcp| fc 2,|| fc 3 : || + k 2,± fc 3,_L 

= -ffo(fe,||,X2;oo)# (fc3,|| ) X3;fcc) + i k 2 ,\\ k 3i \\L (k 2i \\,X2;oo)H (k 3j \\,X3;kc) 



1 - 3 

+ 2 fc2 'll fc 3,||^o(fc2,|| I X2;oo)Lo(/s 3 ,||,X3;fcc) + = A52.ii fe 3 ,|| ^ofe.H , X25 oo)L (k 3 ,\\ , xr, kc) 



1 



- yfc 3 || J?o(fc 2 ,||,X2;oo)L (A:3 i ,|,X3;fcc) - ^fc 2) || L o(k 2 ,\\ , X2; oo)J?o(A5 3 ,|| , X3] kc)- 
The expression for the window function acting on the Lyman-a flux is in this case 



W a (-k 2>ll -k u )W a (k 2i[l ) 



1 - e 



^2, II 



E 



1 - e 
(-2)" (k 



2 1 \ / ft- 



k L 



-4- g fc 3.|]/ fe " (g 2fc2 ||/fe 2 _g 



y(-2)" p2,|iy/fc 3 ,|| 
k J \ k 



(A32) 



(A33) 



where we have recast the window function in a combination that is suitable for furthering the calculation. Notice 
in fact that the first and second term in the sum differ only by the presence of k^ or k in the denominators of the 
exponentials. Furthermore, the terms in square brackets are functions of fc2.11 only. We then define the coefficients 



dkw 



g 2fc|| fk g fc|| /fc 



/ m (fc||,X,oo), 
/ m (fc||,x,oo). 



(A34) 
(A35) 



A point worth making is that the second expression can be obtained from the first one with the substitution kL — ► k 
in the denominators but not in the expression for k, hence the necessity of two separate definitions. Considering then 
the following generic term, it is possible to show that 

dk 2 dk 3lVlQ ~ n _ . \tj r { ,.u rr /^^Ax 



E 



2tt 2tt " 
(-2)" 1 - 



k 2 kl fi(k 2 , X2; oo)gj(k 3 , \3] k c )W a {-k 2 - k 3 ) W a {k 2 ) e 
k (P + ,) g (, +m ) {AXiX . fccj +-)( X2) fcL ) + fc(p+ 9 )^+™)(A x , X ; fc c , k)f^ +m \x2, k)] , (A36) 



which then leads directly to 



(^2,3 = 2 J] 



^St\k L )Ht\A X ; k c , k L ) + \kll^ +1 \k L ) < i+1) (A X ; k c ,k L ) 



1 



2 &(m+l). 



(k L ) 4 m+1) (A X ; fee, fci) + ^kl L^'(k L ) L^'(A X ; k c , k L ) 



4 fO+2), 



-(m+2). 



^^ m) (fc L )4 m+2) (A X ;fc c ,fcL) - ^Lt +2) (k L )Ht\Ax;kc,k L ) + (k L - fc) 



(A37) 



Note that in the above expression while fm are always real, the f m can be real or imaginary depending on whether 
n is even or odd. However, the fact that fm^ is zero whenever the upper index is odd guarantees that (S 2 5)2,3 is always 
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real valued. Also, notice that while the coefficients /m (xq;&l) are decreasing with to, the coefficients fm (xQ>k) 
are actually increasing with to. However, the to! factor present in the denominator more than compensates for these 
increasing coefficients and allows to truncate the series in an actual calculation. 

Finally, it is worth pointing out that the case without cutoff on the long wavelength mode is recovered from the above 
expression simply by setting ki = and then noticing that in this case k = and that therefore the corresponding 
terms appearing in Eqs. (|A29| lA37|) disappear. 
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